An improved determination of the width of the top quark 
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We present an improved determination of the total width of the top quark, Ft, using 5.4 fb-i of 
integrated luminosity collected by the DO Collaboration at the Tevatron pp Collider. The total width 
Ft is extracted from the partial decay width r(t — ^ Wb) and the branching fraction B{t — >■ Wb). 
r{t — >■ Wb) is obtained from the f-channel single top quark production cross section and B{t — ^ Wb) 
is measured in ti events. For a top mass of 172.5 GeV, the resulting width is Ft = 2.00lL°oil GeV. 
This translates to a top-quark lifetime of rt = (3.29lo gg) x 10"'^^ s. We also extract an improved 
direct hmit on the CKM matrix element 0.81 < |V"tb| < 1 at 95% C.L. and a limit of \Vtb>\ < 0.59 
for a high mass fourth generation bottom quark assuming unitarity of the fourth generation quark 
mixing matrix. 

PACS numbers: 14.65. Ha. 14.65. Jk, 12.15.Hh 



INTRODUCTION 

The top quark is the heaviest known elementary parti- 
cle and completes the quark sector of the standard model 
(SM). It differs from the other quarks not only by its 
much larger mass, but also by its lifetime that is ex- 
pected to be shorter than the QCD scale typical of the 
formation of hadronic bound states Within the SM, 
the top quark decays almost exclusively into a W boson 
and a b quark. The total decay width Tt is therefore 
expected to be dominated by the partial decay width 
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r(t — s- Wb). Neglecting higher order electroweak correc- 
tions and terms of order /m^ , a ^ , and (as / tt)M^ / , 
the partial widthprcdicted by the SM at ncxt-to-leading 
order (NLO) is 

where mt (mf,) is the mass of the top (bottom) quark, Gp 
(as) the Fermi (strong interaction) coupling constant, 
Mw the mass of the W boson, and Vtb the strength of 
the left-handed Wtb couphng. Setting as{Mz) = 0.118, 
Gf = 1.16637 X 10~5 GeV-^ Mw = 80.399 GeV, 
l^tbl = 1 HI, and assuming mt = 172.5 GeV, we ob- 
tain r(t Wb)gyi = 1.33 GeV. A deviation from the 
theoretical prediction would indicate the presence of be- 
yond SM (BSM) physics, including those involving BSM 
decays of the top quark to final states that escape detec- 
tion. Examples of such BSM scenarios are the anomalous 
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duction cross section as 

r{t-)- Wb) = ct(^- channel) 



(T(i— channel) 



SM 



The total decay width Ft can be written in terms 
of the partial decay width and the branching fraction 
B{t Wb) as 



FIG. 1: Tree-level Feynman diagrams for (a) s and (b) t- 
channel single top quark production. 



Wtb couplings hadronically decaying charged Higgs 
bosons y or a fourth generation b' quark. 

The electroweak single top quark production at the 
Tevatron proceeds mainly through the exchange of a vir- 
tual W boson accompanied at tree level by a quark in 
the s channel la i or by both a b and a light quark in 
the t channel [y, lZ|. A third channel, tW, in which the 
top quark is produced in association with a W boson, is 
not considered in this analysis because the expected pro- 
duction cross section at the Tevatron is small @. Fig- 
ure [1] shows the tree-level Feynman diagrams for s- and 
i-channel production Q. In this Letter, we determine 
F(i — Wb) from a measurement of the t-channel single 
top quark production cross section, making use of the 
fact that the process involves a Wtb vertex and is thus 
proportional to F(t — Wb). The i-channel was chosen as 
it has the highest production cross section at the Teva- 
tron Q and because BSM contributions may have differ- 
ent effects on the s- and i-channel cross sections. Here 
we do not assume that the s-channel production rate is 
as predicted by the SM. 

The first direct upper bound on F^ [l^l was set by 
the CDF Collaboration from an analysis of the invari- 
ant mass distribution of candidate events using 1 fb^^ 
of integrated luminosity. The first indirect determina- 
tion of Tt was obtained by the DO Collaboration by 
combining the measurement of the single top i-channel 
cross section using 2.3 fb~^ of integrated luminosity 
and the branching fraction B(t — !■ Wb) determined from 
a sample of ti events in 0.9 fb^^ of integrated luminos- 
ity [13[- This method assumes the W ^ tb coupling in 
single top quark production is the same as in top quark 
decay. 

In this Letter, we apply the method of [llj in a new 
indirect determination of Tt that is based on two prior 
DO measurements, both performed using 5.4 fb~^ of in- 
tegrated luminosity: the single top quark t-channel cross 
section a{pp tqb+X) = 2.90±0.59 (stat+syst) pb [H 
and the ratio R = B(t ^ Wb)/B{t Wq) = 0.90 ± 
0.04 [l^l, where q can be a d, s or 6 quark. 

The partial decay width F(t — >■ Wb) = Fp can be ex- 
pressed in terms of the t-channel single top quark pro- 



Ft 



(3) 



B{t Wb) 

Combining Eqs. [3]and[21 the total decay width becomes 



cr(t-channcl) T{t Wb) 



SM 



B{t -> Wb) cr(i- channel) SM ' 



(4) 



from which it is possible to derive the lifetime of the top 
quark as Tt = fi/Tt- 

We can also use the measured value of Fp to probe 
the Wtb interaction and directly determine the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix [l6j 
element \Vtb\- A direct determination of |Vtb|, without 
assuming unitarity of the CKM matrix or three genera- 
tions of quarks, is possible through the measurement of 
the total single top quark production cross section [l7| . 
However, this method assumes that the top quark decays 
exclusively to Wb, and assumes the relative production 
rate of s and t-channel single top production is as pre- 
dicted by the SM. These two assumptions are removed 
when we combine the branching fraction measurement 
(which allows for t — )• Wd and t — >■ Ws decays) and the 
single top cross section measured from the i-channel, in- 
dependently of any assumption on the s-channel rate or 
on the ratio of s- to i-channel production cross sections. 



ANALYSIS METHOD 

This analysis relies on two prior DO measurements, 
the single top t-channel cross section |14| and the ra- 
tio of the top quark branching fraction [l5|. Both are 
based on 5.4 fb~^ of integrated luminosity. The latter 
is performed by distinguishing between the standard de- 
cay mode of the top quark, tt — >■ W~^bW~b, (indicated by 
bb), and decay modes that include light quarks {qi = d,s), 
tt ^ W+bW-qi {bqi) and ti W+qiW'qi (qiqi). The 
analysis relies on a sample of ti events in which one W 
boson decays into a quark and an antiquark and the 
other into an electron or muon and a neutrino, or events 
in which both W bosons decay into liy. In both cases, 
we accept events in which the W boson decayed to a t 
lepton that subsequently decayed into an electron or a 
muon. We use a neural network 6-tagging algorithm 
to identify jets that originate from the hadronization of 
long-lived b hadrons (5-tagged jet) and distinguish be- 
tween the bb, bqi and qiqi final states in ti decay. 
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The t-channel cross section measurement uses events 
containing an isolated electron or muon, missing trans- 
verse energy and at least two jets. Background is sup- 
pressed by requiring that one or two of the jets is identi- 
fied as a 6-jet. The main background contributions arise 
from W bosons produced in association with jets and 
from tt pairs. We further improve the discrimination be- 
tween signal and background by em ploy ing multivariate 
analysis techniques as described in [19|. We use a dis- 
criminant trained to separate the ^-channel signal from 
the backgrounds in 6 independent analysis channels, de- 
fined according to jet multiplicity (2, 3 or 4), and number 
of 6-taggcd jets (1 or 2) [Ti |. 

Based on the t-channel output discriminant distribu- 
tion, wc define a binned likelihood 



M 



L(Did)=n' 



A! 



(5) 



where cr^ is the s-channel cross section times B{t ^ Wb), 
SLt and are arrays containing the product of the ac- 
ceptance and the integrated luminosity in each bin for t 
and s processes, respectively, and b is an array contain- 
ing the mean count of expected background events. The 
term c^p t} is given by 



B{t M^6)(T(t-channel)sM 

r(t ^ wb)sM 



or by 



Ct 



(7) 



B'^it W6)cr(t-channel) 



SM 



Tit ^ Wb)sM 



(8) 



where D and d are arrays containing the observed and 
mean expected count for all M bins from the six different 
analysis channels. The mean expected count can be writ- 
ten in terms of the partial (Tp) or total (Fj) top quark 
width as 

d(r{p^t},cr^,at,as,b) = c^p^t^Tip^t^at + a'^a^ +b (6) 



when measuring the partial or total top quark decay 
width, respectively. The extra B{t ^ Wb) term with re- 
spect to Eqs.[2]and|4]is needed to remove the assumption 
of S(t — )• Wb) ~ 1 used when generating the single top 
quark and tt samples. We then form a Bayesian proba- 
bility density for the partial or total width by integrating 
the expression 



p(r{p,t}|D) = J L{B\T ^p^t} , a'^, at, as,h)Tr{T ^p^t})T^{a'jTr{at)n{as)TT{h)dasdatda.sdh, 



(9) 



where 7r(*) represent our prior knowledge of the parame- 
ters a'g, at , as and b . The normalization constant Af 
ensures that / p(T^p t}\T>)dT^p ty = 1. The integration is 
performed assuming a positive and uniform probability 
density for and for as- The other priors quantify 

our knowledge of the systematic uncertainties for the val- 
ues of at, as and b. Each independent systematic source 
is modeled with a Gaussian of mean zero and width set to 
one standard deviation of the corresponding uncertainty. 



SYSTEMATIC UNCERTAINTIES 

The main systematic uncertainties affect the t-channcl 
output discriminant as well as the measured branching 
fraction i3(i — J- Wb), and are summarized in Table HI 
Common systematics that affect both the discriminant 
and B{t — ^ Wb) are taken as 100% correlated. 

The terms included in the uncertainty calculation are: 

• Uncertainty on jet flavor identification involving b, 
c and light-fiavor jet tagging rates and the calorime- 



Sources 



Size 



Uncertainties on B{t — ^ Wb) 

b-jet identification 4.0 

tt cross section 2.1 

Integrated luminosity 1.6 

Statistical uncertainty 2.3 

Statistical correlation 4.2 
Uncertainties on a{t — channel) 

6-jet identification 9.3 

tt cross section 3.1 

Integrated luminosity 5.1 

W-I-Jets normalization 8.1 

Jet energy resolution 11.6 

Jet energy scale 6.8 

Monte Carlo statistics 6.7 



TABLE I: Sources of statistical and main systematic uncer- 
tainties relative to the measured value for t-channel cross sec- 
tion and branching fraction that affects the determination of 
the partial/total decay width. We list the the most impor- 
tant uncertainties for the branching fraction and t-channel 
cross section measurements, respectively. 
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tcr response to &-jets. 

• Uncertainties on the modeling of the process used 
in the B(t — )■ Wb) measurement, including uncer- 
tainties from parton distribution functions (PDFs) , 
different event generators and hadronization mod- 
els. These uncertainties are correlated with the tt 
background yield uncertainty in the i-channel dis- 
criminant. 

• Uncertainty on the integrated luminosity. 

• Uncertainties on the normalization of the VF-l-jets 
heavy-flavor contribution. 

• Uncertainty in the jet energy scale and resolution. 

• The statistical uncertainty of the B{t Wb) mea- 
surement is added as systematic uncertainty. The 
statistical uncertainty related to the t-channel cross 
section is by construction included in the top quark 
width posterior distributions. 

• A systematic uncertainty is added to account for 
the statistical correlation due to a small overlap 
between events selected to build the t-channel dis- 
criminant and those used in the B{t — > Wb) mea- 
surement. 

• Signal and background yield uncertainty because of 
the amount of Monte Carlo events used to construct 
the background and signal discriminant. 

Other terms included in the calculation of i-channel 
discriminant but not mentioned in the table are: 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
r(t ^ Wb) [GeV] 

FIG. 2: [color online] Probability density for the expected and 
measured partial width Tp. The shaded areas represent one 
standard deviation centered on the most probable value. 
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+0.39 
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GeV 



0.5 1 1.5 



.= 2.00 GeV 



2.5 3 3.5 4 4.5 
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FIG. 3: [color online] Probability density for the expected 
and measured total width Ft. The shaded areas represent 
one standard deviation around the most probable value. 



• Uncertainties on modeling the single top quark 
signal, including initial- and final-state radiation, 
scale uncertainties and PDFs. 

• Detector simulation uncertainty arising from the 
modeling of particle identification in the simulated 
samples. 

• Uncertainties arising from the modeling of the dif- 
ferent background sources that are obtained using 
data-driven methods. 



quark width is found to be 

F, = 2.00t"oil GeV, (11) 
which can be expressed as a top quark lifetime of 

(12) 



This also translates in an upper limit to the top quark 
lifetime of n < 4.88 x IQ-^s s at 95% C.L.. 



RESULT 

The expected and observed probability densities for 
the partial width Fp are shown in Fig. [2] The most prob- 
able value for the partial width is defined by the peak of 
the probability density function and corresponds to 



Fp = 1.87t^:l^ GeV. 



(10) 



The expected and observed probability densities for 
the total width Ft are shown in Fig. |31 The total top 



TOP QUARK MASS DEPENDENCE 

The measured branching fraction and f-channel pro- 
duction cross section, as well as T{t — > Wb)sM, depend 
on the top quark mass mj. To study this dependency, 
we repeat the analysis using simulated tt and single top 
samples that were generated at different values of mt in 
the range 170 to 175 GeV. The value of F(t ^ Wb)sM 
is recalculated depending on mt- Given that the depen- 
dence from mt is small, the value and uncertainties for 
B{t Wb) corresponding to nit = 172.5 GeV are used 
in all cases. 
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mt (GeV) 


170 


172.5 


175 


r{t Wb)sM (GeV) 
r(t ^ W6) (GeV) 
Tt (GeV) 


1.26 


1.33 

1 C7+0.44 
^•UU_o 43 


1.40 

1 yf.+0.40 
i.yD_0.43 



TABLE II: Observed partial and total top quark decay width 
as a function of the top quark mass. We provide also the 
values for SM top quark partial widths used in the analysis. 



Tabic im summarizes the variation of the partial and to- 
tal top quark decay width as a function of TOj . The table 
also lists the values of r{t — > Wb)sM used in the analy- 
sis. The variation of the decay width with rrit follows the 
non-monotonic variation observed for the t-channel cross 
section 



The effect of the mass dependency can be quantified 
by interpolating the observed Tt in function of top mass 
from table |TT] to the current Tevatron combination mt = 
173.2 ± 0.9 GeV Adding a mass uncertainty of 



STt = max \rt{mt) 

mte[172.5,175] 



rt(173.2)| « 0.07 GcV 



in quadrature to the symmetrized interpolated error for 
irit = 173.2 GeV results in a value for the total width of 
Tt = 2.03 ± 0.46 GeV, and a value for top quark lifetime 
of Tt 3.24l[] J^-10"2^ s. A lower limit on the total decay 
width Ft > 1.37 GeV at 95% C.L. can be estimated by 
assuming that the posterior density for Tt approximates 
a Gaussian distribution. This translates into an upper 
limit on the top quark lifetime of Tt < 4.82 • 10^^^ s at 
95% C.L. Therefore we conclude that the effect of the 
mass uncertainty is small with respect to the observed 
uncertainty obtained assuming a top quark mass of mt = 
172.5 GcV. 



MEASUREMENT OF |Vtb| 

Wc construct a posterior probability density for \Vtb\ 
by setting 

C{p.t}'^{p,t} = \Vtb\'^B{t-^ Wb)a{t~cha.nncl)sM,\Vt,\=i 

(13) 

in Eq.El A lower limit of |Vt6| > 0.81 at the 95% C.L. is 
obtained by restricting the prior for \Vtb\'^ to be uniform 
for < |Vtbp < 1, as illustrated in Fig. S) A systematic 
uncertainty on the theoretical prediction for the t-channel 
cross section was included in the result. 

We apply the same procedure to constrain the strength 
of the coupling of a fourth generation b' quark to the 
top quark |Vt6'|. For this measurement we assume that 
mb' > mt — mw, a small probability density for the b' 
quark to exist in protons and antiprotons, and unitarity 
of the four-generation quark-mixing matrix with \ Vtb\'^ + 
\Vtb'\^ = 1, and \Vtd\,\Vt,\ 1. Using the limit on 
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|V I > 0.81 at 95% C.L. 
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FIG. 4: [color online] Distributions of the posterior density for 
jVtbp. The shaded (darked shaded) band indicates regions of 
68% (95%) C.L. 



|Vtf,| and the condition \Vtb'\'^ = 1 — |Vt6p, wc obtain 
|ytf,'| < 0.59 at the 95% C.L. 



SUMMARY 

Wc have presented an improved determination of the 
width of the top quark using the Baycsian techniques 
previously used to measure the single top quark pro- 
duction cross section and an improved measurement of 
the branching fraction B{t Wb). The method assumes 
that the coupling leading to t-channcl single top quark 
production is identical to the coupling in the top quark 
decay. We have determined the top quark width as Ft = 
2.00toil GeV for mt = 172.5 GeV, which corresponds to 
a top quark hfetime of n = (3.29^°;^°) x lO'^'^ s. These 
are the most precise determinations of width and lifetime 
to date. In addition, we set a lower limit of |Vth| > 0.81 at 
the 95% C.L. without assuming that the top quark decays 
exclusively to Wb and with no assumption on the s- and 
t-channcl relative production rate. Wc also set a limit on 
the strengths of the coupling for a fourth-generation b' 
quark to the top quark of |Vtfc'| < 0.59 at 95% C.L. 
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